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The consequences of hadro-chemical freezeout for the subsequent hadron gas evolution in central 
heavy-ion collisions at RHIC and LHC energies are discussed with special emphasis on effects due 
to antibaryons. Contrary to naive expectations, their individual conservation, as implied by experi- 
mental data, has significant impact on the chemical off-equilibrium composition of hadronic matter 
at collider energies. This may reflect on a variety of observables including source sizes and dilepton 
spectra. 
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The systematic study of hadron production in central 
collisions of heavy nuclei at (ultra-) relativistic energies 
has revealed strong evidence for a high degree of equili- 
bration achieved in the course of these reactions. In par- 
ticular, it has been shown that, over a large range of colli- 
sion energies (from SIS to SPS) [|l|-§, the measured final 
state hadron abundances can be well described by es- 
sentially two parameters, the temperature T and baryon 
chemical potential fiB (which bear little sensitivity on 
centrality 0|). Importantly, the inferred collision-energy 
dependence reinforces confidence in the notion that one is 
actually probing different regions of the QCD phase dia- 
gram in these experiments. The analyses are in line with 
the naive expectation that with increasing CM-energy 
y/s of the colliding nuclei, the central rapidity regions 
are increasingly governed by interactions between \ow-x 
partons of the incoming nucleons, thus leading to a de- 
crease in the net baryon content. This is reflected in a 
systematic decrease of the baryon chemical potential with 
•y/s, accompanied by an increase in temperature. In fact, 
for the top SPS energies {^/s = 17 AGeV) and for the 
first RHIC data at ^/s = 130 AGeV, the such extracted 
(/i, T) values are very close to the expected phase bound- 
ary to the Quark-Gluon Plasma (QGP). Under RHIC 
conditions this implies a copious production of baryon- 
antibaryon pairs. Whether those are solely of thermal 
origin, or what exactly the underlying microscopic (par- 
tonic) production mechanism is, remains a matter of de- 
bate at present. It is clear, however, that, as with any 
hadronic observable, the impact of the later hadronic 
stages has to be well understood before any firm con- 
clusions can be drawn. 

At SPS and RHIC/LHC energies the total hadron den- 
sities implied by the hadro-chemical analyses are not 
small, and one cannot expect the system to decouple at 
this point. It has been realized, however, that hadronic 
cross sections for elastic scattering, especially when in- 
volving intermediate resonances {e.g., tttt p, ttN <-!■ 
A, ttK ^ K*), are typically much larger than those 
for inelastic (number-changing) reactions. The associ- 
ated timescales for thermal equilibration are thus much 
smaller than for the chemical one. At SPS energies and 
above, thermal equilibrium in the expanding hadron gas 
can be maintained for about 10 fm/c. This is much be- 
low typical chemical relaxation times and has led to the 



picture of a chemical freezeout at {^ichiTch) (as inferred 
from hadron ratios) with Tch — 160-180 MeV, followed by 
a thermal decoupling at significantly lower temperatures, 
Tth — 100-120 MeV. The existence of this intermediate 
hadronic phase is well supported as it importantly fig- 
ures into both hadronic H] and electromagnetic observ- 
ables §. 

An immediate consequence of this picture is that all 
hadron species which are not subject to strong decays, 
have to be effectively conserved in number subsequent 
to chemical freezeout. In statistical mechanics language, 
this can be described by the build-up of (effective) chem- 
ical potentials, /z^, fiK, Pr], etc., in the hadronic evo- 
lution towards thermal freezeout. Corresponding ther- 
modynamic trajectories at SPS energies have been con- 
structed ]20|-|2^. In this context, antibaryons play a spe- 
cial role. Although their production from AGS to RHIC 
energies follows chemical-freezeout systematics, this can- 
not be expected a priori, since their annihilation cross 
section on baryons is large {e.g., app^nn — 50 mb with 
an average n=5-6 at the typical thermal energies). Under 
SPS conditions, this issue has been resolved in ref. [p5| , 
where it is shown that the backward reaction of multi- 
pion fusion, in connection with the rather large oversatu- 
ration of the pion phase space towards thermal freeze- 
out (entailing /ij/* ~ 70 MeV), is capable of approx- 
imately maintaining the chemical-freezeout number of 
antibaryons. On the other hand, with a p/p ratio of 
5-7% [|l^,ll[ and a final n/B ratio (after strong decays) 
of around 5, the detailed treatment of antibaryons does 
not exert a large infiuence on the bulk thermodynam- 
ics of the hadronic evolution. The main point of this 
note is to show that this is no longer true at RHIC ener- 
gies. Based on the experimental fact that also at RHIC 
the chemical freezeout abundances survive the hadron 
gas phase, we will argue that the explicit conservation 
of antibaryon number significantly affects the composi- 
tion (and possibly evolution) of the hadronic expansion 
at collider energies. In particular, we will construct an 
explicit isentropic trajectory in the (J-b-T plane including 
effective number-conservation constraints. 

Before we proceed with a more quantitative assess- 
ment, let us further motivate our objective by the follow- 
ing estimate: RHIC data at y/s = 130 AGeV for central 
Au-Au collisions p2|-p4[ imply a midrapidity density of 
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charged pions plus kaons of about dNch/dy ~ 600, so that 
dN^r+K /dy ~ 900. On the other hand, the total baryon 
number has been estimated to be dNs/dy ~ 10 [|l5|Jl6[ 
so that dNgj^g/dy ~ 165 (using B/B ~ 0.65 flJlll)- 
Since a (anti-) baryon carries about twice the entropy 
than that of pion, it follows that the former account for 
at least 25% of the total entropy at midrapidity. 

Along the lines of standard hadro-chemical analyses 
we base the following calculations for the thermody- 
namic state variables on an ideal hadron gasQ includ- 
ing a set of resonance states comprised of 37 meson and 
37 baryon species (up to masses tom = 1.7 GeV and 
ruB = 2 GeV , respectively). We will focus on conditions 
representative for midrapidities in central Au+Au colli- 
sions at ^/s = 200 AGeV (very similar results emerge for 
a net-baryon free environment appropriate for LHC). For 
definiteness, we employ a specific entropy (i.e., entropy 
per net baryon) of S/N^^*^ = 250 which, together with 
a chemical freezeout temperature of T^h = 180 MeV, re- 
sults in a baryon chemical potential of /j,^ = 24 MeV. 
This agrees well with (moderate) extrapolations of very 
similar thermal model analyses |^,^ based on RHIC 
data at y/s = 130 AGeV. By definition, all meson- 
chemical potentials are zero at Tch, and we also neglect 
(small) isospin (^7) and strangeness (^s) chemical poten- 
tials. As usual, the effective pion number, which is sub- 
ject to conservation, is defined by including all strongly 
decaying resonances with lifetimes shorter than the du- 
ration of the interacting hadronic phase, i.e.. 



7V;?^/ = F7^B^7V«ft(T,^,) 



(1) 



and likewise for if , 77 and 77'. In eq. (|i|), Vfb denotes the 
fireball 3- volume (which is related to the centrality of the 
collision but docs not affect thermodynamic quantities) 
and Qi the number density of hadron i. The effective 
pion-number Tvi'' of a given resonance is determined by 
the number of pions in its decay modes including the 

branching ratios, e.g., M''' = 2, N^'^^ = M^''' = 1, 
^^^(1520)) ^ Q ^ ^ _|_ Q^^^ ^ 2. The same weighting 

applies to the corresponding effective chemical poten- 
tial of each resonance, e.g., /ip = 2^^^, /^A = ^J■N + /^tt, 
tJ-K' = fJ-n- + fix, MAf(i520) ^ ^J■N + l-45^7r, which mcaus 
that the strong decays and formations are in relative 
chemical equilibrium. For a few states this distinction 
is not very sharp, most notably uj and (j) mesons. Their 
free lifetimes are well above hadronic fireball lifetimes, 
but in-medium effects are expected to alter this situa- 
tion 1 28 - 33 1 . In the former case, some of the w's and 
0's decay after chemical freezeout without being regen- 
erated, i.e., their number decreases; in the latter case. 



*for a discussion of the influence of hadronic in-medium ef- 
fects, see, e.g., refs. Iislhsl]. 



their abundances stay (at least for a while) in chem- 
ical equilibrium which, in the case of large pion- and 
kaon-chemical potentials, could even increase their num- 
bers. In our calculation we employ an intermediate sce- 
nario keeping the individual uj and (j) numbers approxi- 
mately constant (which is numerically close to the chem- 
ical equilibrium assumption). An important point now 
is the treatment of antibaryons as first pointed out in 
ref. [Q. We also conserve their number by introducing 
another effective chemical potential /il/'^; this implies. 
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FIG. 1. Isentropic thermodynamic trajectories at RHIC 
(left) and SPS (right) energies, starting from the empiri- 
cally determined chemical- freezeout points |l|J^j5|j26| . Full 
lines include antibaryon-number conservation, dashed lines 
do not. (Also note that the trajectories converge towards the 
zero-temperature axis at /ijv = mjv.) 

A thermodynamic trajectory in the fJ-N-T plane is now 
constructed starting from chemical freezeout under the 
additional assumption of total entropy conservation (as 
applied in ideal fluid dynamics). Together with net 
baryon number conservation (which is, of course, ex- 
act) this amounts to keeping the specific entropy per net 
baryon fixed, S/N^"^ = s/g'^''* = 250 for RHIC with 
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d^k 



[±/(c^,)ln/(c^,) 



(2) 



(27r)3 

+ (lT/(^0)ln(lT/(c^»)] 

/^•(a.b„T)-/^-(A7b.,T) (3) 



d^k 



(27r)3 



denoting entropy- and net baryon number-density, re- 
spectively (in eq. (|^), upper (lower) signs refer to 
fermions (bosons), and / are the pertinent thermal dis- 
tribution functions). The resulting trajectory is shown 
by the full line in Fig. |l|. For comparison, we also dis- 
play a trajectory without antibaryon-number conserva- 
tion (dashed line). The difference is rather dramatic. 
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indicating a strongly increased baryon chemical poten- 
tial towards the regions of expected thermal freezeout at 
around Tth — 120 MeV. This is induced by an approxi- 
mately linear increase of the antinuclcon chemical poten- 
tial, which is displayed in Fig. ^ along with nucleon, pion, 
kaon, and eta chemical potentials. At chemical freezeout 
the sum of baryon and antibaryon density is close to nu- 
clear saturation density po — 0.16 fm^'^, Qb+b — 1-lpo- 
On the other hand, at T = 120 MeV, Qb+b - ^-"^Bq 
for the solid-line trajectory, but Qb+b — 0.02po for the 
dashed-line trajectory. 

One can readily verify that the solid-line trajectory in 
Fig. |l| is consistent with the measured p/p ratio, which 
is approximately given by exp[— A/iat/T] with A// at — 
fiM — I^N- It is also important to note that in the yf^^^ = 
case, the pion-chemical potentials remain close to zero. 
Thus, B conservation at collider energies is at the origin 
of large meson chemical potentials in the late hadronic 
stageaj. The reason can be traced back to the substantial 
amount of entropy stored in BB excitations. At a given 
temperature, the smaller available entropy per pion at 
fixed pion number can only be realized by a finite , as 
also noted in rcf. Q. To illustrate this fact, we display 
in Fig. ^ the entropy per pion in a ir-p gas. One finds 
that for T = 120 MeV, S/N^ is reduced by about 20% 
when raising /ijr from zero to its actual value of ~80 MeV 
corresponding to the solid-line RHIC-trajectory of Fig. |^. 
This is in line with the rough estimate of the entropy 
carried by BB pairs at RHIC above. 
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FIG. 2. Temperature dependence of nucleon-, pion-, kaon-, 
eta- and antinucleon-chemical potentials for an isentropic 
hadronic fireball expansion at RHIC energy (200 AGeV). 
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is opposite to SPS conditions, where, following the 
arguments of ref. pion-oversaturation is responsible for 
regenerating antibaryons. 
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FIG. 3. Entropy per pion in a n-p gas at fixed temperature 
as a function of pion-chemical potential (using pp = 2^^ and 
counting p's as two pions). 

It furthermore follows that at given temperature the 
fireball volume Vp b of the system with finite p^^^ is sub- 
stantially reduced, i.e., the expanding hadron gas cools 
faster. E.g., at T = 120 MeV the ratio of the 3-D firebaU 
volumes for the = a-^d p'^^^ > solution for the 
trajectories amounts to a factor of ~3. Whether this will 
translate into smaller source size radii as extracted via 
HBT correlation measurement (and thus contribute to an 
explanation of the current 'HBT-puzzle' [p5[) remains to 
be investigated within a (hydro-) dynamical simulation. 
Other observable consequences of the increased densities 
in the hadronic phases at RHIC concern low-mass dilep- 
ton emission, where in-medium effects are rather sensitive 
to (the sum of) anti-/baryon densities |^,^. Accordingly 
improved estimates have been presented in ref. ||3^ . 

In our analysis we did not allude to underlying mecha- 
nisms that could provide a constant antibaryon abun- 
dance in the hadronic phase at RHIC. As mentioned 
above, at SPS energies this issue could be resolved p^ ] 
within a thermal framework by accounting for multi-pion 
fusion reactions in the presence of a large Pt^. Although 
it is conceivable that this mechanism is also operative at 
RHIC, the situation is more involved due to the feedback 
of the B abundance on Pt^. This problem is not easily 
addressed within transport or cascade simulations either 
due to the difficulties of treating the backward reaction of 
multi-pion fusion (see ref. for recent progress). Cor- 
responding trajectories that have been extracted from, 
e.g., UrQMD transport calculations are therefore re- 
sembling the dashed lines in Fig. 

In summary, based on the notion of hadro-chemical 
freezeout, we have constructed (isentropic) thermody- 
namic trajectories for the hadronic phase in ultra- 
relativistic heavy-ion collisions. Consistency with the 



3 



observed particle ratios has been enforced by effective 
conservation laws, implemented via finite chemical po- 
tentials for stable hadrons [w.r.t. strong interactions). 
The main emphasis was on consequences of conserved 
antibaryon-number. At SPS energies and below, these 
are small, but at collider energies (RHIC and LHC) we 
have found appreciable modifications in the composition 
of the expanding hadronic system. In particular, the en- 
tropy carried by BB pairs is not small and has been iden- 
tified as the origin of large meson-chemical potentials, 
which in turn leads to smaller and thus denser systems 
at temperatures prior to thermal freezeout. This could 
have important consequences for hadronic and electro- 
magnetic observables, such as source radii and medium 
effects in low-mass dilepton production. The former are 
rather sensitive to the detailed expansion dynamics of a 
heavy-ion reaction, which we did not address here. For 
the same reason, the question where on the given trajec- 
tories thermal freezeout in the hadronic phase occurs has 
not been answered, but needs to be investigated in hy- 
drodynamical and/or microscopic transport simulations. 

This work was supported by the U.S. Department of 
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